To be or not to be: orb, the fusome and oocyte specification in Drosophila.

In the fruit fly Drosophila melanogaster, two cells in a cyst of 16 interconnected cells have the potential to become the oocyte, but only one of these will assume an oocyte fate as the cysts transition through regions 2a and 2b of the germarium. The mechanism of specification depends on a polarized microtubule network, a dynein dependent Egl:BicD mRNA cargo complex, a special membranous structure called the fusome and its associated proteins, and the translational regulator orb. In this work, we have investigated the role of orb and the fusome in oocyte specification. We show here that specification is a stepwise process. Initially, orb mRNAs accumulate in the two pro-oocytes in close association with the fusome. This association is accompanied by the activation of the orb autoregulatory loop, generating high levels of Orb. Subsequently, orb mRNAs become enriched in only one of the pro-oocytes, the presumptive oocyte, and this is followed, with a delay, by Orb localization to the oocyte. We find that fusome association of orb mRNAs is essential for oocyte specification in the germarium, is mediated by the orb 3' UTR, and requires Orb protein. We also show that the microtubule minus end binding protein Patronin functions downstream of orb in oocyte specification. Finally, in contrast to a previously proposed model for oocyte selection, we find that the choice of which pro-oocyte becomes the oocyte does not seem to be predetermined by the amount of fusome material in these two cells, but instead depends upon a competition for orb gene products.

A dynamical model of growth and maturation in Drosophila.

The decision to stop growing and mature into an adult is a critical point in development that determines adult body size, impacting multiple aspects of an adult's biology. In many animals, growth cessation is a consequence of hormone release that appears to be tied to the attainment of a particular body size or condition. Nevertheless, the size-sensing mechanism animals use to initiate hormone synthesis is poorly understood. Here, we develop a simple mathematical model of growth cessation in Drosophila melanogaster, which is ostensibly triggered by the attainment of a critical weight (CW) early in the last instar. Attainment of CW is correlated with the synthesis of the steroid hormone ecdysone, which causes a larva to stop growing, pupate, and metamorphose into the adult form. Our model suggests that, contrary to expectation, the size-sensing mechanism that initiates metamorphosis occurs before the larva reaches CW; that is, the critical-weight phenomenon is a downstream consequence of an earlier size-dependent developmental decision, not a decision point itself. Further, this size-sensing mechanism does not require a direct assessment of body size but emerges from the interactions between body size, ecdysone, and nutritional signaling. Because many aspects of our model are evolutionarily conserved among all animals, the model may provide a general framework for understanding how animals commit to maturing from their juvenile to adult form.

A multifunctional Wnt regulator underlies the evolution of rodent stripe patterns.

Animal pigment patterns are excellent models to elucidate mechanisms of biological organization. Although theoretical simulations, such as Turing reaction-diffusion systems, recapitulate many animal patterns, they are insufficient to account for those showing a high degree of spatial organization and reproducibility. Here, we study the coat of the African striped mouse (Rhabdomys pumilio) to uncover how periodic stripes form. Combining transcriptomics, mathematical modelling and mouse transgenics, we show that the Wnt modulator Sfrp2 regulates the distribution of hair follicles and establishes an embryonic prepattern that foreshadows pigment stripes. Moreover, by developing in vivo gene editing in striped mice, we find that Sfrp2 knockout is sufficient to alter the stripe pattern. Strikingly, mutants exhibited changes in pigmentation, revealing that Sfrp2 also regulates hair colour. Lastly, through evolutionary analyses, we find that striped mice have evolved lineage-specific changes in regulatory elements surrounding Sfrp2, many of which may be implicated in modulating the expression of this gene. Altogether, our results show that a single factor controls coat pattern formation by acting both as an orienting signalling mechanism and a modulator of pigmentation. More broadly, our work provides insights into how spatial patterns are established in developing embryos and the mechanisms by which phenotypic novelty originates.

A model of replicating coupled oscillators generates naturally occurring cell networks.

When a founder cell and its progeny divide with incomplete cytokinesis, a network forms in which each intercellular bridge corresponds to a past mitotic event. Such networks are required for gamete production in many animals, and different species have evolved diverse final network topologies. Although mechanisms regulating network assembly have been identified in particular organisms, we lack a quantitative framework to understand network assembly and inter-species variability. Motivated by cell networks responsible for oocyte production in invertebrates, where the final topology is typically invariant within each species, we devised a mathematical model for generating cell networks, in which each node is an oscillator and, after a full cycle, the node produces a daughter to which it remains connected. These cell cycle oscillations are transient and coupled via diffusion over the edges of the network. By variation of three biologically motivated parameters, our model generates nearly all such networks currently reported across invertebrates. Furthermore, small parameter variations can rationalize cases of intra-species variation. Because cell networks outside of the ovary often form less deterministically, we propose model generalizations to account for sources of stochasticity.

Phase plane dynamics of ERK phosphorylation.

The extracellular signal-regulated kinase (ERK) controls multiple critical processes in the cell and is deregulated in human cancers, congenital abnormalities, immune diseases, and neurodevelopmental syndromes. Catalytic activity of ERK requires dual phosphorylation by an upstream kinase, in a mechanism that can be described by two sequential Michaelis-Menten steps. The estimation of individual reaction rate constants from kinetic data in the full mechanism has proved challenging. Here, we present an analytically tractable approach to parameter estimation that is based on the phase plane representation of ERK activation and yields two combinations of six reaction rate constants in the detailed mechanism. These combinations correspond to the ratio of the specificities of two consecutive phosphorylations and the probability that monophosphorylated substrate does not dissociate from the enzyme before the second phosphorylation. The presented approach offers a language for comparing the effects of mutations that disrupt ERK activation and function in vivo. As an illustration, we use phase plane representation to analyze dual phosphorylation under heterozygous conditions, when two enzyme variants compete for the same substrate.

Dynamics of an incoherent feedforward loop drive ERK-dependent pattern formation in the early Drosophila embryo.

Positional information in development often manifests as stripes of gene expression, but how stripes form remains incompletely understood. Here, we use optogenetics and live-cell biosensors to investigate the posterior brachyenteron (byn) stripe in early Drosophila embryos. This stripe depends on interpretation of an upstream ERK activity gradient and the expression of two target genes, tailless (tll) and huckebein (hkb), that exert antagonistic control over byn. We find that high or low doses of ERK signaling produce transient or sustained byn expression, respectively. Although tll transcription is always rapidly induced, hkb converts graded ERK inputs into a variable time delay. Nuclei thus interpret ERK amplitude through the relative timing of tll and hkb transcription. Antagonistic regulatory paths acting on different timescales are hallmarks of an incoherent feedforward loop, which is sufficient to explain byn dynamics and adds temporal complexity to the steady-state model of byn stripe formation. We further show that 'blurring' of an all-or-none stimulus through intracellular diffusion non-locally produces a byn stripe. Overall, we provide a blueprint for using optogenetics to dissect developmental signal interpretation in space and time.

Evolutionary history of MEK1 illuminates the nature of deleterious mutations.

Mutations in signal transduction pathways lead to various diseases including cancers. MEK1 kinase, encoded by the human MAP2K1 gene, is one of the central components of the MAPK pathway and more than a hundred somatic mutations in the MAP2K1 gene were identified in various tumors. Germline mutations deregulating MEK1 also lead to congenital abnormalities, such as the cardiofaciocutaneous syndrome and arteriovenous malformation. Evaluating variants associated with a disease is a challenge, and computational genomic approaches aid in this process. Establishing evolutionary history of a gene improves computational prediction of disease-causing mutations; however, the evolutionary history of MEK1 is not well understood. Here, by revealing a precise evolutionary history of MEK1, we construct a well-defined dataset of MEK1 metazoan orthologs, which provides sufficient depth to distinguish between conserved and variable amino acid positions. We matched known and predicted disease-causing and benign mutations to evolutionary changes observed in corresponding amino acid positions and found that all known and many suspected disease-causing mutations are evolutionarily intolerable. We selected several variants that cannot be unambiguously assessed by automated prediction tools but that are confidently identified as "damaging" by our approach, for experimental validation in Drosophila. In all cases, evolutionary intolerant variants caused increased mortality and severe defects in fruit fly embryos confirming their damaging nature. We anticipate that our analysis will serve as a blueprint to help evaluate known and novel missense variants in MEK1 and that our approach will contribute to improving automated tools for disease-associated variant interpretation.

Maintaining symmetry during body axis elongation.

Bilateral symmetry defines much of the animal kingdom and is crucial for numerous functions of bilaterian organisms. Genetic approaches have discovered highly conserved patterning networks that establish bilateral symmetry in early embryos,1 but how this symmetry is maintained throughout subsequent morphogenetic events remains largely unknown.2 Here we show that the terminal patterning system-which relies on Ras/ERK signaling through activation of the Torso receptor by its ligand Trunk3-is critical for preserving bilateral symmetry during Drosophila body axis elongation, a process driven by cell rearrangements in the two identical lateral regions of the embryo and specified by the dorsal-ventral and anterior-posterior patterning systems.4 We demonstrate that fluctuating asymmetries in this rapid convergent-extension process are attenuated in normal embryos over time, possibly through noise-dissipating forces from the posterior midgut invagination and movement. However, when Torso signaling is attenuated via mutation of Trunk or RNAi directed against downstream Ras/ERK pathway components, body axis elongation results in a characteristic corkscrew phenotype,5 which reflects dramatic reorganization of global tissue flow and is incompatible with viability. Our results reveal a new function downstream of the Drosophila terminal patterning system in potentially active control of bilateral symmetry and should motivate systematic search for similar symmetry-preserving regulatory mechanisms in other bilaterians.

Self-organized intracellular twisters.

Life in complex systems, such as cities and organisms, comes to a standstill when global coordination of mass, energy, and information flows is disrupted. Global coordination is no less important in single cells, especially in large oocytes and newly formed embryos, which commonly use fast fluid flows for dynamic reorganization of their cytoplasm. Here, we combine theory, computing, and imaging to investigate such flows in the Drosophila oocyte, where streaming has been proposed to spontaneously arise from hydrodynamic interactions among cortically anchored microtubules loaded with cargo-carrying molecular motors. We use a fast, accurate, and scalable numerical approach to investigate fluid-structure interactions of 1000s of flexible fibers and demonstrate the robust emergence and evolution of cell-spanning vortices, or twisters. Dominated by a rigid body rotation and secondary toroidal components, these flows are likely involved in rapid mixing and transport of ooplasmic components.

Self-organized intracellular twisters.

Life in complex systems, such as cities and organisms, comes to a standstill when global coordination of mass, energy, and information flows is disrupted. Global coordination is no less important in single cells, especially in large oocytes and newly formed embryos, which commonly use fast fluid flows for dynamic reorganization of their cytoplasm. Here, we combine theory, computing, and imaging to investigate such flows in the Drosophila oocyte, where streaming has been proposed to spontaneously arise from hydrodynamic interactions among cortically anchored microtubules loaded with cargo-carrying molecular motors. We use a fast, accurate, and scalable numerical approach to investigate fluid-structure interactions of 1000s of flexible fibers and demonstrate the robust emergence and evolution of cell-spanning vortices, or twisters. Dominated by a rigid body rotation and secondary toroidal components, these flows are likely involved in rapid mixing and transport of ooplasmic components.

Dynamics of an incoherent feedforward loop drive ERK-dependent pattern formation in the early Drosophila embryo.

Positional information in developing tissues often takes the form of stripes of gene expression that mark the boundaries of a particular cell type or morphogenetic process. How stripes form is still in many cases poorly understood. Here we use optogenetics and live-cell biosensors to investigate one such pattern: the posterior stripe of brachyenteron (byn) expression in the early Drosophila embryo. This byn stripe depends on interpretation of an upstream signal - a gradient of ERK kinase activity - and the expression of two target genes tailless (tll) and huckebein (hkb) that exert antagonistic control over byn . We find that high or low doses of ERK signaling produce either transient or sustained byn expression, respectively. These ERK stimuli also regulate tll and hkb expression with distinct dynamics: tll transcription is rapidly induced under both low and high stimuli, whereas hkb transcription converts graded ERK inputs into an output switch with a variable time delay. Antagonistic regulatory paths acting on different timescales are hallmarks of an incoherent feedforward loop architecture, which is sufficient to explain transient or sustained byn dynamics and adds temporal complexity to the steady-state model of byn stripe formation. We further show that an all-or-none stimulus can be 'blurred' through intracellular diffusion to non-locally produce a stripe of byn gene expression. Overall, our study provides a blueprint for using optogenetic inputs to dissect developmental signal interpretation in space and time.

A novel ground truth dataset enables robust 3D nuclear instance segmentation in early mouse embryos.

For investigations into fate specification and cell rearrangements in live images of preimplantation embryos, automated and accurate 3D instance segmentation of nuclei is invaluable; however, the performance of segmentation methods is limited by the images' low signal-to-noise ratio and high voxel anisotropy and the nuclei's dense packing and variable shapes. Supervised machine learning approaches have the potential to radically improve segmentation accuracy but are hampered by a lack of fully annotated 3D data. In this work, we first establish a novel mouse line expressing near-infrared nuclear reporter H2B-miRFP720. H2B-miRFP720 is the longest wavelength nuclear reporter in mice and can be imaged simultaneously with other reporters with minimal overlap. We then generate a dataset, which we call BlastoSPIM, of 3D microscopy images of H2B-miRFP720-expressing embryos with ground truth for nuclear instance segmentation. Using BlastoSPIM, we benchmark the performance of five convolutional neural networks and identify Stardist-3D as the most accurate instance segmentation method across preimplantation development. Stardist-3D, trained on BlastoSPIM, performs robustly up to the end of preimplantation development (> 100 nuclei) and enables studies of fate patterning in the late blastocyst. We, then, demonstrate BlastoSPIM's usefulness as pre-train data for related problems. BlastoSPIM and its corresponding Stardist-3D models are available at: blastospim.flatironinstitute.org.

Fusome topology and inheritance during insect gametogenesis.

From insects to mammals, oocytes and sperm develop within germline cysts comprising cells connected by intercellular bridges (ICBs). In numerous insects, formation of the cyst is accompanied by growth of the fusome-a membranous organelle that permeates the cyst. Fusome composition and function are best understood in Drosophila melanogaster: during oogenesis, the fusome dictates cyst topology and size and facilitates oocyte selection, while during spermatogenesis, the fusome synchronizes the cyst's response to DNA damage. Despite its distinct and sex-specific roles during insect gametogenesis, elucidating fusome growth and inheritance in females and its structure and connectivity in males has remained challenging. Here, we take advantage of advances in three-dimensional (3D) confocal microscopy and computational image processing tools to reconstruct the topology, growth, and distribution of the fusome in both sexes. In females, our experimental findings inform a theoretical model for fusome assembly and inheritance and suggest that oocyte selection proceeds through an 'equivalency with a bias' mechanism. In males, we find that cell divisions can deviate from the maximally branched pattern observed in females, leading to greater topological variability. Our work consolidates existing disjointed experimental observations and contributes a readily generalizable computational approach for quantitative studies of gametogenesis within and across species.

Stochastic phenotypes in RAS-dependent developmental diseases.

Germline mutations upregulating RAS signaling are associated with multiple developmental disorders. A hallmark of these conditions is that the same mutation may present vastly different phenotypes in different individuals, even in monozygotic twins. Here, we demonstrate how the origins of such largely unexplained phenotypic variations may be dissected using highly controlled studies in Drosophila that have been gene edited to carry activating variants of MEK, a core enzyme in the RAS pathway. This allowed us to measure the small but consistent increase in signaling output of such alleles in vivo. The fraction of mutation carriers reaching adulthood was strongly reduced, but most surviving animals had normal RAS-dependent structures. We rationalize these results using a stochastic signaling model and support it by quantifying cell fate specification errors in bilaterally symmetric larval trachea, a RAS-dependent structure that allows us to isolate the effects of mutations from potential contributions of genetic modifiers and environmental differences. We propose that the small increase in signaling output shifts the distribution of phenotypes into a regime, where stochastic variation causes defects in some individuals, but not in others. Our findings shed light on phenotypic heterogeneity of developmental disorders caused by deregulated RAS signaling and offer a framework for investigating causal effects of other pathogenic alleles and mild mutations in general.

Quantitative models for building and growing fated small cell networks.

Small cell clusters exhibit numerous phenomena typically associated with complex systems, such as division of labour and programmed cell death. A conserved class of such clusters occurs during oogenesis in the form of germline cysts that give rise to oocytes. Germline cysts form through cell divisions with incomplete cytokinesis, leaving cells intimately connected through intercellular bridges that facilitate cyst generation, cell fate determination and collective growth dynamics. Using the well-characterized Drosophila melanogaster female germline cyst as a foundation, we present mathematical models rooted in the dynamics of cell cycle proteins and their interactions to explain the generation of germline cell lineage trees (CLTs) and highlight the diversity of observed CLT sizes and topologies across species. We analyse competing models of symmetry breaking in CLTs to rationalize the observed dynamics and robustness of oocyte fate specification, and highlight remaining gaps in knowledge. We also explore how CLT topology affects cell cycle dynamics and synchronization and highlight mechanisms of intercellular coupling that underlie the observed collective growth patterns during oogenesis. Throughout, we point to similarities across organisms that warrant further investigation and comment on the extent to which experimental and theoretical findings made in model systems extend to other species.

Mechanics of stabilized intercellular bridges.

Numerous engineered and natural systems form through reinforcement and stabilization of a deformed configuration that was generated by a transient force. An important class of such structures arises during gametogenesis, when a dividing cell undergoes incomplete cytokinesis, giving rise to daughter cells that remain connected through a stabilized intercellular bridge (ICB). ICBs can form through arrest of the contractile cytokinetic furrow and its subsequent stabilization. Despite knowledge of the molecular components, the mechanics underlying robust ICB assembly and the interplay between ring contractility and stiffening are poorly understood. Here, we report joint experimental and theoretical work that explores the physics underlying robust ICB assembly. We develop a continuum mechanics model that reveals the minimal requirements for the formation of stable ICBs, and validate the model's equilibrium predictions through a tabletop experimental analog. With insight into the equilibrium states, we turn to the dynamics: we demonstrate that contractility and stiffening are in dynamic competition and that the time intervals of their action must overlap to ensure assembly of ICBs of biologically observed proportions. Our results highlight a mechanism in which deformation and remodeling are tightly coordinated-one that is applicable to several mechanics-based applications and is a common theme in biological systems spanning several length scales.

Putting in the Erk: Growth factor signaling and mesoderm morphogenesis.

It has long been known that FGF signaling contributes to mesoderm formation, a germ layer found in triploblasts that is composed of highly migratory cells that give rise to muscles and to the skeletal structures of vertebrates. FGF signaling activates several pathways in the developing mesoderm, including transient activation of the Erk pathway, which triggers mesodermal fate specification through the induction of the gene brachyury and activates morphogenetic programs that allow mesodermal cells to position themselves in the embryo. In this review, we discuss what is known about the generation and interpretation of transient Erk signaling in mesodermal tissues across species. We focus specifically on mechanisms that translate the level and duration of Erk signaling into cell fate and cell movement instructions and discuss strategies for further interrogating the role that Erk signaling dynamics play in mesodermal gastrulation and morphogenesis.

Dynamics of Drosophila endoderm specification.

During early Drosophila embryogenesis, a network of gene regulatory interactions orchestrates terminal patterning, playing a critical role in the subsequent formation of the gut. We utilized CRISPR gene editing at endogenous loci to create live reporters of transcription and light-sheet microscopy to monitor the individual components of the posterior gut patterning network across 90 min prior to gastrulation. We developed a computational approach for fusing imaging datasets of the individual components into a common multivariable trajectory. Data fusion revealed low intrinsic dimensionality of posterior patterning and cell fate specification in wild-type embryos. The simple structure that we uncovered allowed us to construct a model of interactions within the posterior patterning regulatory network and make testable predictions about its dynamics at the protein level. The presented data fusion strategy is a step toward establishing a unified framework that would explore how stochastic spatiotemporal signals give rise to highly reproducible morphogenetic outcomes.

Deconstructing gastrulation at single-cell resolution.

Gastrulation movements in all animal embryos start with regulated deformations of patterned epithelial sheets, which are driven by cell divisions, cell shape changes, and cell intercalations. Each of these behaviors has been associated with distinct aspects of gastrulation1-4 and has been a subject of intense research using genetic, cell biological, and more recently, biophysical approaches.5-14 Most of these studies, however, focus either on cellular processes driving gastrulation or on large-scale tissue deformations.15-23 Recent advances in microscopy and image processing create a unique opportunity for integrating these complementary viewpoints.24-28 Here, we take a step toward bridging these complementary strategies and deconstruct the early stages of gastrulation in the entire Drosophila embryo. Our approach relies on an integrated computational framework for cell segmentation and tracking and on efficient algorithms for event detection. The detected events are then mapped back onto the blastoderm shell, providing an intuitive visual means to examine complex cellular activity patterns within the context of their initial anatomic domains. By analyzing these maps, we identified that the loss of nearly half of surface cells to invaginations is compensated primarily by transient mitotic rounding. In addition, by analyzing mapped cell intercalation events, we derived direct quantitative relations between intercalation frequency and the rate of axis elongation. This work is setting the stage for systems-level dissection of a pivotal step in animal development.

ERK signaling dissolves ERF repression condensates in living embryos.

Phase separation underlies the organization of the nucleus, including the biogenesis of nucleoli and the packaging of heterochromatin. Here we explore the regulation of transcription factor condensates involved in gene repression by ERK signaling in gastrulating embryos of a simple proto-vertebrate (Ciona). ERK signaling induces nuclear export of the transcriptional repressor Ets-2 repressive factor (ERF), which has been linked to various human developmental disorders. Using high-resolution imaging, we show that ERF is localized within discrete nuclear condensates that dissolve upon ERK activation. Interestingly, we observe dynamic pulses of assembly and dissociation during interphase, providing visualization of a nuclear phase separation process regulated by cell signaling. We discuss the implications of these observations for producing sharp on/off switches in gene activity and suppressing noise in cell-cell signaling events.